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Abstract. We derive a general expression for the Hankel determinants of a 
Dirichlet series F{s) and derive the asymptotic behavior for the special case 
that F(s) is the Riemann zeta function. In this case the Hanliel determinant 
is a discrete analogue of the Selberg integral and can be viewed as a matrix in- 
tegral with discrete measure. We briefly comment on its relation to Plancherel 
measures. 
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1. Introduction 
In this paper we will consider the Hankel determinants 
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and various generalizations of them. These determinants go very rapidly to zero as 
the dimension of the matrix becomes large e.g. 



4.9 X 10-16684 
4.3 X 10-16871 



1 2 7 1 16 1 

'2n+ 1 ^ (2n+ 1)2 ~ 3 (2n + l)3 ^ T (2n + 1)^ 



41 

¥(2n 



The ratios have (experimentally) a surprisingly simple asymptotic expansion: 

Hl^'lml^UC] ^ 1 1,21 6 1 56 1 
i?(°)[C]HP[C] 2n (2n)2^3(2n)3 5 (2n)4 ^ 45 (2n)5 ^ ' ' ' ' 

This recursion can be solved to yield explicit expressions for H^\(]. In fact more 
detailed numerical experiments by D. Zagier confirmed our findings with the result 
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240 (2n)2 ~ 7257600 (2n)4 ~ 1741824000 (2n)6 
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with the interesting observation that 

« 0.351466738331 

^9/8 

Wc do not know how to prove the full asymptotic expansions but will describe 
a method that lets one at least understand the weaker asymptotic form 

logH(°)[C] ~ logif(i)[C] ~ -n2(log(2n) - 

We will also discuss interesting relations with a continuous version (Selberg inte- 
gral) and with the Plancherel measure of the symmetric group. We start with the 
presentation of some results on the Hankel determinants of Dirichlet series. 

2. Hankel determinants of Dirichlet series 
Definition. Let F{s) be a Dirichlet series with coefficients /(n) i.e. 

^' fin) 



, n" 



For n integer, n > 0, r integer, r > we define the Hankel determinant Hn\F]: 

ifW[i^]=det(i^(i + i + r)),<. .<„ 
We will also use the notation Hn[F] = H^\f]. Our first result is: 
Theorem 2.1. Hn^ [F] is given by 



mi,m2,---Tnn — l i—1 ^ i<j 



Proof. We prove the case r = first. Using the definition of the determinant 

H^^^F] = J2(-^)''F{l + 7r{l))F{2 + n{2))...F{n + n{n)) 
7ree„ 

/(mi)/(m2)---/(m„) 



E E (-^y 



l+7r(l) 2+7r(2) n+n{n) 



mi...m„>l7reS„ ™i "^'> ■ ■ ■ m„ 



mi...mn>i i ^ " ^gg^ m2 ...m„ 
^ /(mi)/(m2)---/(mn) -pr / 1 1_ 

mi...mn>l i ^ i<j ^ J 

where is the symmetric group and (—1)'^ is the sign of the permutation tt. In 
the last line we have used the Vandermonde determinant 

X- til! = I nf--- 
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Interchanging two summation variables say and rrij with i ^ j the Vander- 
monde determinant changes sign. Summing over all permutations of {l,2,...n} 
and dividing by the number of permutations we obtain: 



^(0) \p^^L /(toi)/(to2) ■■■/(m„) {-lY \\ ( — 

mi...m„>l 7re6„ ^\^> 7r(n) v 
^ J_ f (mi) f (1712) ■■ ■/(?•»„) (-1)'' TT / J- 

n! ^ mim2---rn„ ^ ™'^(i)™'^(2) _7r(n) 11 ^ ^ 

■mi...7nn>l ttGC 
^ j_ /(toi)/(to2) ■■■ fjnin) -r-r 

nl ^-^ mirni ■ ■ ■ vn^ 

mi...m„>l 1 z (t 

1 /(™i)/(™2) •••/(m„) -pp 
— r > 7 I I {'mi — m,) 



For r > we replace f{n) f{n)n ^ which gives 

ifW[F] = l V /(m,)/(m2).--/K0 -pr 2 
" ^ ^ n! ^ (toiTO2---?7t-«)^"+'^ 11^ * ^' 

mi...m,i>l i<j 

which completes the proof. □ 



Corollary 2.2. If f{n) is multiplicative it follows trivially from Theorem 2.1 



n! (mim2 ■ - . mn) ^ 

mi,m2,...m„ — 1 z<j 

which is again a Dirichlet series. Note that 

and vanishes if the number of prime factors of m is less than n—\ and is divisible by 
iXYi^ii'i ~ 1)0^ since each summand Wi^j {mi — mj)^ is divisible by (Jir=i(* ~ -'^)')- 
For a proof see j2] . 

Lemma 2.3. If f{n) > for all integers n > then 

Proof. The smallest m contributing in the sum is m = 1 • 2 • • • n = n!. For rrii = 
{1, 2, ... n} we have 



n n 



mi-m2---m„=m i<j i=\ j=i+l \i=l / 



SO there is at least one term greater than zero in the sum in Corollary |2.2[ The 
other terms in the sum are all greater or equal to zero so the sum is positive which 
proves the lemma. □ 
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Definition 2.4. For each integer n > and integer m > we define the function 

K{m) by 

hn{m) ^ — ^ (TOi - TOj)^ . 

For n = 2 this function can be expressed in terms of arithmetic functions: 

h2{m) ~ (72(1^) — md(m) 

where a'2{m) is the divisor function <T2{m) — J2d\m'^'^ ^^'^ d{m) = J2d\m gives the 
number of divisors of m. We are now in the position to state our second theorem. 

Theorem 2.5. For a Dirichlet series F{s) with coefficients multiplicative coeffi- 
cients / the Hankel determinant Hn[F] is given by a Dirichlet series with coefficients 
h„{m)f{'in): 

hn{m)f{m) 



Proof. Inserting the definition of hn{m) in Corollary 2.2 gives the right hand ex- 
pression. □ 

Example 2.6. Our first example, already discussed in the introductions, is f{n) = 
1, so that F{s) — C(s) where ({s) is the Riemann zeta function. Using Theorem 
Iwe obtain: 



(2.1) = l y TT (m, - m,)' . 

■mi...inn>l i<j 

In particular Hn\C] > for any n. Specializing to e.g. n — 2 and r = we obtain 
from Theorem 12.11 



mi,m2>l 
1 ^ 



2 V mimi m^m^ ra^nd 

mi,m2>l \ J- ^ J- J- ^ 

C(2)C(4)-C(3)' 

C(2) C(3) 
C(3) C(4) 

The Hankel determinant can also be expressed as a linear combination of multiple 
zeta values 

Hn[(\ = ^ (-1) ^ l+7r(l) 2+7r(2) n+7r(n) 

Tree,, m„>...m2>mi>l "^1 "^2 ' ' ' ™n 

= E (-irC(l + 4l),2-f7r(2),...,n-H^(n)). 

TreSn 

Example 2.7. The second nontrivial example is 

^{m) , 1 ^ TJ , n2 



det 



771 — 1 \ in-i-ni2---mn—in i<^j 

where ii{n) is the Mobius function. 



+ i) J l<^.J<r * 
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Remark. Multiple Dirichlet series have been investigated for some time for a recent 
work see e.g. |3] and |3]. 

3. A PROBABILISTIC MODEL FOR THE ASYMPTOTIC BEHAVIOR OF Hn[C] 

In this section we determine the behavior of i?n[C] as rt ^ oo. The basic idea 
is to find the dominant contribution to the sum. We note that all contributions 
are positive and that the Vandermonde determinant is only nonzero if the nii are 
pairwise different. We can reorder them so that to„ > . . . > 1712 > mi. There are 
precisely n\ contributions with toi,TO2 • • ■ "^n in the unordered sum. We write 

(3.1) Hn[C] = V exp($(mi,TO2 . . .m„)) 



m ^ > . . . > m2 > m 1 > 1 



with 



(3.2) <[>(mi, m2 . . . m„) = — log(mi) + 2 ^ ^ log (m^ — m^) . 

i—l i—1 

Finding the largest ^{1711,1712 .. .nin) is a discrete combinatorial optimization 
problem. We can also view $ as an energy functional of a one-dimensional Coulomb 
gas problem on a lattice where a large (— 2n) attractive charge is placed at zero 
and the charges placed at positions integer repel each other with a logarith- 
mic potential. Unlike in the standard one-dimensional Coulomb gas problem, the 
charges cannot have a distance smaller than one. Let us consider the case n ^ 1. 
Building up the configuration by adding a charge at one by one for the first few 
mi, 1712 ■ ■ ■ the first term is dominant so that $ is optimized by placing the first 
charge at mi = 1, the second at m2 = 2 and so on. Adding more charges slowly 
builds up the second term which produces a repulsive potential which makes it more 
favorable to place charges at > i. From this analogy we expect the density of 
the mi to be one from i = 1 up to some imax < n and then to decay to zero as 
i ^ n. 

We define the distribution function for each configuration {mi, 7712, . . . m„} 



(3.3) p{x) = 5{nx — mj) 

i=l 

with 5{x) being the delta distribution. The integral over p{x) is 

/■oo 

(3.4) / p{x)dx = 1. 







by definition. We call a distribution which obeys Eq. |3.4| normalized. The func- 
tional <i>(mi, TO2 . . . m„) can be expressed as a functional of the distribution function 
p{x): 



£ 



^[p] = — n (2 / p{x)log{nx)dx — // log \n{x ~ y)\p{x)p{y)dxdy 

-n^ ( 21og(n) -I- 2 / p{x)log{x)dx — log{n) + // p{x) p{y) log \x — y\dxdy 
-n"^ log(?T.) — n"^ [ 2 / p{x)log{x)dx — // p{x) p{y) log \x — y\ dxdy 



(3.5) 
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We seek a continuous test function p{x) which maximizes the functional Eq. 3.5 
subject to the constraint that < p{x) < 1 for all x and p{x) = for a; < and 
Eq. |3.4[ We assume that such a continuous functions exists. The arguments for its 



existence can probably be made much more rigorous by using Poissonization, see 



e.g. jT] 13 . The constraint Eq. 3.4 can be taken care of by introducing a Lagrange 
multiplier A € M and finding an extremum of <i>[/ci] — A p[x]dx. The problem of 
finding the extremum of ^[p] is thus reduced to the integral equation 

nOO 

(3.6) 2 log(a;) - 2 / p{y) log\x - y\ ^ A. 

This integral equation only applies when p{x) can actually be varied, i.e. for p{x) ^ 
1 and p{x) ^ 0. Differentiating with respect to x eliminates the Lagrange multiplier 
and finally leads to 

(3.7) 1 = f ^dy 

X Jo x-y 

where ^ denotes the principal value integral. Assume p{x) < 1 for all a; € IR then 
the integral equation has to be fulfilled for all x G [0,cx)). The integral equation 
has as only solution p{x) = S{x) where 6{x) is the Dirac distribution which diverges 
for a; ^ and does not fulfill the constraint. Therefore the constraint p{a) < 1 has 
to be sharp for some positive real number a. We will use an Ansatz for p{x) and 
verify that it obeys the integral equation and the constraint. 



Theorem 3.1. Define the function p{x) for x E [0,oo): 



(0 < X < i) 



Then p{x) is the continuous solution of the stationary condition Eq. 3.7 for a; > 1/2 



is normalized and fulfills the condition 1 > p(x) > for all positive x G 
Proof. The derivative of p is: 



(3.8) p'{x) 



ix<l/2) 
Let e G K with e > 0. Define ^^(a:) as 



(3.9) A,{x) = i / log I (a: - yj' + | p'{y)dy. 







where log(2) denotes the principal branch of the logarithm (0 < arg(z) < tt). Using 
and substituting y = {s'^ + l)/2 we have 



3.8 



^ ^ 4 p I log I (s^ + (1 - 2a;)) + - log(2) 



log 



(s2 + (i_2x))' + 



ds + log(2) 



HANKEL DETERMINANTS OF DIRICHLET SERIES 



7 



First consider the case x > 1/2. Using contour integration above the real axis 
closing the contour in the upper half plane gives 



2m — 
as 



1 log ((.s2 + 1 - 2xf + 
vr 



-2i 



1 



4s(s2 + 1 - 2a;) 
(1 + s2 _ 2a;) V e2 {s + i)^ 

T - log 



-2 



+ log(2) 

-i 

log((l + g2-2a;)' + e2) 
(s + i)3 



log 



4 



Taking the limit e ^ we have 



lim AJx) 



log\x\ 



Using integration by parts on Eq. |3.9| we have for all e > 



A, 



:{x) = ^ log ({x ~ yf + e^^ p{x) 



x—oo 
x=l/2 



(x - y) 



In the limit e ^ we have 



lim A^(x) — — log(x) 



1/2 (a; — y) + 

p{y) 



p{v)dy 



-dy 



where f denotes the principal value integral. Combining the two expressions for 
A^{x) we obtain 



1 

X 



p{x) 



lo x~y 

for a; > 1/2 which is the first statement of the theorem. 
The normalization integral is 1 since 



p{x)dx = ^ + 
1 

= 2 + 



1/2 



arctan 



1 



1 



ii/2 r 



V2a;- 1 
arctan(y) — 



1 + r 



^/2x-l 

2i 
dy 



dx 



1 1 arctan(y) 



y=oo 



1 



y=i/2 A/2 2/(1 + 



dy 



The function p{x) is continuous at 1/2 since for e e M with ^ > e > 

and p {-^ — £) = 1- The function p{x) is continuous and monotonically decreasing 
with increasing x since p'{x) < for x G (1/2, oo). The maximum of p{x), x G 
[0, oo) is one and the infimum is since for x G [0, oo) with a; ^ 1 

p(^)^^a;-3/2 + ofa;-^/2) > 0. 
which completes the proof. □ 
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To determine ^[p] we need to evaluate the integral p{y) log \x — y\. For x > 
1/2 we can use the stationary condition Eq. 3.6 however we have to determine A 
first. 



Proposition 3.2. The function p{x) x G (0,oo) of Theorem 3.1 has the property 



log \x - y\ p{y)dy = < 



log (2) - 1 {x = 0) 

-VI - 2a; - log(2) + 2{x - 1) log {l + ^J I - 2x) + x\og{2x) (0 < a; < 1/2) 
[log(x) (1/2 < x) 



Proof. Consider the case a; > 1/2 first. From theorem |3.1| ve have 



p{y) log \x -y\dy = log(x) + -A. 



To prove A = it is sufficient to prove it for one a; > 1/2. We choose x — 1/2. 
Using the definition of p{x) we have 



p{y) log 



1 



dy = 



1/2 



log 





log 
log 
log 
log 



1 

2 " 
1\ 1 



dy+ I (log(2y-l)-log(2))p(2/)dy 

1/2 



2J 2 

l\ 1 
2 

1 
2 
1 
2 



1/2 



\og{2y~ l)p{y)dy 



1/2 



{{2y-l)\og{2y-l)^{2y~l))p'{y)dy 



1 2 



oo 9 2 



2s^ log(s) - 



2 TT io (s2 + lY 



ds 



where we have substituted s = l/\/2x — 1 .The integral over s can be done by 
contour integration. Finally we have 



p{y) log 



This implies A = which proves the case for x > 0. 

Next consider the case < a; < 1/2 . We split the integral in two parts and 
integrate by parts 



p{y) log \x - y\ 



\dy 



1/2 /.oo 

log \x-y\dy+ I p{y) log (y - x) i 

J 1/2 

-X + X \og{x) + / {{y ~ x) log {y - x) - {y - x)) dy 

Jl/2 
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The second integral can be done by substituting s = \/2y — 1 . The second term 
in the last line is 



P'(y)(2/ - ^) (log (y-x)- 1) dy 



1/2 



2 ("^ 

71" io 



log (s2 + 1 - 2.t) - log(2) - 1 



(s2 + l)^ 



(s-^ + 1 - 2x)ds 



(log(2) + 1) (x - 1) + 



log (s^ + 1-2, 



(s2 + 1) 



2 (s^ + 1 - 2a;)ds 



= (log(2) + l)(a;-l) + 21og(l + Vr^^) + 

where the last integral was evaluated using contour integration (see e.g. |5], 4.295, 
integral 7). Collecting all terms we obtain the a; < 1/2 case of the proposition. 
Finally consider the limit x ^ with a: > 0. We have to leading order in x 

-Vl - 2x-log(2)+2(x-l) log (1 + Vl - 2a;)+a;log(2a;) = -l+log(2)+0(xlog(a;)). 

In the limit a; ^ 0, a; > we find 

/■oo 

lim / p{y)\o%\x - y\dy = -l+log(2) 

2:^0+ Jo 

which completes the proof. □ 



1 + VI - 2a: 



We are now in the position to evaluate (t)[p\. Using Proposition 3.2 all integrals 
of Eq. 3.1 can be reduced to elementary integrals with the result 

<i> [p] = — log(n) — ^2 J p{x) log{x) dx — JJ p{x) p{y) log \x — y\ dxdy^ 

= -n^ log(n) - (2 (log(2) - 1) - f log(2) 

= -n^ ^log(2n) - 
From this we conclude that the dominant contribution to Hn{C) is 

log(iJ„(C))«-n2 (^log(2,i)-0 
which agrees with numerical findings. 



4. Relation to the Selberg integral 

We next discuss the relation of Hn\C) to the Selberg integral (Selberg's exten- 
sion of the beta integral [TIJ, for a detailed explanation see [1] chapter 8) which 
plays a central role in random matrix theory (see |T0], chapter 17) and is given by 

5„(a, /3, 7) = tr' ■ ■ ■ ' (1 - tif^' • • • (1 - tnf^' n - ^^)'^ '^^i • • • 

^0 Jo 
Substituting ti l/nii we write 



Sn{a,(3,'-f) = 
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For a— l + r, P~l and 7 = 1 we find 

(4.1p„ (r + 1, 1) = / •••/ „ TT {nii - nijf dmi ■ ■ ■ dnin 

Ji Ji (mi---TO„) 

The similarity between Eq. |4.1| and Eq. |2.1| is striking and can be generalized easily. 
It can be seen immediately from the definition that S'„(l, 1, 1) > for each n > 0. 

However as we will prove now for n ^ 1 the Selberg integral ^^(l, 1, 1) is much 
larger than i?„(C) excluding a naive application of Euler-MacLaurin summation 
formula to i?„(C). It is instructive to repeat the saddle point analysis of the previous 
chapter for S'„(l, 1, 1) in the limit n ^ 00. First note that 



^n(l,l,l) 

with the density 

PS 



/OO POO 
■ J exp{^[ps])dmi ■ ■ ■ drrir, 



1 " 

(x) = - S{x - rrii). 



n 

i=l 



Note that we did not rescale x. The density is normalized to one 

ps{x)dx — 1. 



/ 

J — ( 



Proposition 4.1. The asymptotic density pg maximizing <i>[p] is 

'0 x<\ 



Ps{x) 

X\/X—1 

Proof. The condition that <i>[p] is stationary is 



x> I 



X Ji x-y 

the only difference to Eq. 3.6 is that here no constraint ps{x) < 1 has to be imposed 
since there is no restriction for the difference \mi — mj \ for two integration variables. 
The integral equation can again be solved by standard methods and yields the 
result stated above. Evaluating <i>[ps] we find: 

*[ps] = -n^ [2 J ps{x)log{x)dx- JJ ps{x) ps{y) log \x - y\dxdy^ 

= -n2(41og(2)-21og(2)) 
= -21og(2)n2 

giving S'„(l, 1, 1) « exp(-21og(2)n2 + . . .) so that 5„(1, 1, 1) > i/„(C) > 0. This is 
in fact the correct behavior as we will prove now. □ 

Proposition 4.2. The asymptotic behavior of S'„(l, 1, 1) as n 00 is 

log (^„(1, 1, 1)) - -2 log (2) + (log(27rn) - 1) n + 0(1). 

Proof. The asymptotic behavior can be derived from the exact expression (Theorem 
(8.1.1) in [JQ) 

n 

iog(5„(i,i,i)) = ^(2iog(r(j)) + iog(r(j + i))-iog(r(n + j))) 

= 3 log {G{n +!))+ log (G(n + 2)) - log {G{2n + 1)) . 



HANKEL DETERMINANTS OF DIRICHLET SERIES 



11 



where T{x) is the Euler Gamma function and G{n) — Y[7=o f*-"" n ^ N, n > 
is the Barnes function. Inserting the asymptotic expansion of the Barnes function 



- + - log (27r) z 

A 2 ' 12 



log(z) + 0(l) 



in the previous expression we find 

log(^„(l, 1, 1)) = -21og(2)n2 + (log(27rn) - l)n + 0(1). 
which completes the proof. 

5. Relation to the Plancherel measure 



□ 



We finally remark on the relation of our results to the asymptotic behavior of 
Plancherel measures. The Plancherel measure is defined as 

fdimXY I\l<^<J<n i^^-mjf 

where dim A is the dimension of the representation of the symmetric group S (|A|) 
indexed by the partition A. The sum over all partitions is given by 



E 



1 



(mi ! • • • m„ !)2 

0<mi<...<m„ ^ ' i<j 



An analysis similar to the one above yields the following integral equation for the 
density 

ppiy) 



log(a;) 



-dy. 



lo x-y 

In this case the solution has finite support and is given by 



_ j i arccos (1 - a;) a;e[0,2] 







otherwise 



In this case the constraint pp (0) < 1 does not restrict the solution because pp (0) — 1 
and pp{x) is monotonically decreasing for a; > 0. Below we show for comparison 
the density p{x) (solid line) and pp{x) (dashed line). 




^ X 



Integrating pp{x) over x we recover the famous asymptotic behavior of Plancherel 
measure |12l |9] . 
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